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ABSTRACT 



We present Far Ultraviolet Spectroscopic Explorer (FUSE) observations toward four stars in the 
Large Magellanic Cloud (LMC) superbubble N70 and compare these spectra to those of four com- 
parison targets located in nearby field and diffuse regions. The N70 sight lines show O VI A1032 
absorption that is consistently stronger than the comparison sight lines by ~ 60%. We attribute the 
excess column density (iVovi = iQ^^ os cm~^) to hot gas within N70, potentially the first detection 
of O VI associated with a superbubble. In a survey of 12 LMC sight lines, Howk et al. (2002a) 
concluded that there was no correlation between ISM morphology and iVovi • We present a reanalysis 
of their measurements combined with our own and find a clear difference between the superbubble 
and field samples. The five superbubbles probed to date with FUSE show a consistently higher mean 
TVovi than the 12 non-superbubble sight lines, though both samples show equivalent scatter from halo 
variability. Possible ionization mechanisms for N70 are discussed, and we conclude that the observed 
O VI could be the product of thermal conduction at the interface between the hot, X-ray emitting gas 
inside the superbubble and the cooler, photoionized material making up the shell seen prominently in 
Ha. We calculate the total hydrogen density uh implied by our O VI measurements and find a value 
consistent with expectations. Finally, we discuss emission-line observations of O VI from N70. 
Subject headings: galaxy: halos — ISM: bubbles — magellanic clouds — ISM: structure — supernova rem- 
nants 



1. INTRODUCTION 

Superbubbles (SBs) are some of the most dramatic ex- 
amples of the interaction between hot stars and the sur- 
rounding interstellar medium (ISM). They usually ap- 
pear as limb-brightened rings ~ 100 pc in diameter in 
optical emission lines, often surrounding stellar OB asso- 
ciations. The combined stellar winds of the most massive 
association members excavate a bubble in the surround- 
ing ISM and form a characteristic shell-like structure. 
These bubbles are further heated and energized as the 
most massive association members reach the end of their 
lives and explode as supernovae (SN). The coronal gas 
in the centers of SBs often appears bright in soft X-rays 
while the outer shell is bright in optical emission lines. 
However, little is known observationally about the inter- 
face between the hot interior and the cool shell. Mea- 
suring the hot gas content of SBs will help constrain the 
hot gas production in the ISM, disk-halo feedback mech- 
anisms, as well as refine the models of the life-cycle of 
SBs. 

The advent of the Far Ultraviolet Spectroscopic Ex- 
plorer {FUSE!) satellite has made it possible to ob- 
serve the O VI resonance doublet (AA1031.926, 1037.617), 
which is a sensitive probe of T = 10^ — 10^ K gas. 
Since O^^ has an ionization potential of 114 eV, pho- 
toionized O VI is negligible in most ISM contexts. 
Moderate-velocity shocks, however, can produce O VI 
(fsVinn k ^ 150 km s~^ ^artigan. Ravmond. & Har tmaiml 
1198711 . Gas in the - 10^ K phase can also be 
formed at the interface between hot, coronal gas 
{T > 10^ K) and cooler material through thermal 
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The Magellanic Clouds contain a number of classic 
examples of SBs. However, nearly every sight line to- 
ward the Magellanic Clouds shows O V I absorption at 
both Galactic and Magellanic velocities l |panfort h et al] 
I200I IHoones et aLll2002t IHowk et aI.ll2n02aV For in- 
stance, iHowk et all ()2002aD measured the Large Magel- 
lanic Cloud (LMC) O VI absorption in the spectra of 
twelve carefully-chosen stars and found that A^ovi has 
very little correlation with warm-disk ISM morphology 
as traced by Ha or soft X-ray brightness. Secondly, there 
are variations in Nqyi column density by a factor of four 
or more between sight lines and these variations occur 
over relatively small angular scales. The O VI absorption 
profiles show little relationship to tracers of cooler ions 
and are usually centered at approximately —30 km 
in relation to warm-disk probes such as Fe II. Further- 
more, the O VI profiles are all much broader than can be 
accounted for by simple thermal broadening, implying 
that they sample multiple hot gas systems or turbulent 
absorbers along each line of sight. These results are gen- 
erally consistent with the idea that the observed O VI 
at LMC velocities is due primarily to a halo of hot gas 
surrounding the LMC and expanding toward us up to 
several 10s of km s~^. 

In this paper, we use FUSE observations of a carefully 
selected set of stars within and adjacent to the LMC 
SB N70 l|Henizdll956|) to disentangle the LMC halo and 
SB hot gas signatures. The LMC is a nearly face-on 
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Fig. 1. — N70 is a textbook example of an isolated superbubble, as seen in this Ha+[N II] image. Four stars observed with the FUSE 
MDRS aperture lie clustered in the center of a limb-brightened ring ~105 pc in diameter. Boxes in the upper right show the locations of 
the corresponding serendipitous 30"x30"LWRS aperture positions during the stellar observations. (See text.) 



TABLE 1 

Log of fuse N70 and Comparison Observations 



Target 


RA (J2000) Dec 
h m s ° ' '/ 


N70 
dist. 


Spectral 
Type 


V 

mag 


Fioso" 


FUSE 
ID 


Exp 
(ksec) 


Notes 


LH 114-7 


05 43 12.85 -67 41 16.2 




02Iab 


13.66 


120 


D09811 


13.9 


N70 


Sk-67°250 


05 43 15.48 -67 51 09.6 




09III 


12.68 


150 


D09812 


12.5 


N70 


D 301-1005 


05 43 08.30 -67 50 52.4 




09V 


14.11 


50 


D09814 


3.1 


N70 


D301-NW8 


05 43 15.96 -67 49 51.0 




07V 


14.37 


40 


D09815 


12.2 


N70 


BI272 


05 44 23.18 -67 14 29.3 


37' 


07: II-III:'' 


13.28 


130 


P11729 


10.9 


DEM 310 


Sk-67°266 


05 45 52.00 -67 14 25.0 


39' 


08:Iab: 


12.01 


170 


B09002 


4.9 


DEM 310 


BI 237 


05 36 14.68 -67 39 19.3 


41' 


03V 


13.86 


25 


E51140 


35.6 


diffuse HIT 


CAL 83 


05 43 34.20 -68 22 22.0 


31' 


CV/XRB 


16.20 


1.5 


B01501 


45.8<= 


field 


''Approximate continuum flux at A ^ 1050 A in 


units of (10~ 


-'^ ergs 


cm~^ s 


-1 A-i). 







iMassa et al. 12003) note uncertainty in spectral type; wind lines suggest this star may be hotter and closer 
to the main sequence. 

'^LiFl and LiF2 data were coadded for this sight line for an effective exposure time of ~ 90 ksec. 



disk galaxy (i = 35 ± 6°, Ivan der Marel fc Cioni 



at a well-known distance 51 kpc. iCionr^Tld 



2001) 



200C) 



and free from the line-of-sight confusion and relative dis- 
tance uncertainties present when comparing most Galac- 
tic sight lines. 



We 



2. OBSERVATIONS & DATA ANALYSIS 
chose N70 (D EM 30 1, IHenizd 1195^ 



iDavies. Elhott. fc Meaburnl 119761) as our ideal can 
didate SB to probe for hot gas. It is a beautiful 
ring nebula ^ 105 pc in diameter centered 2.6° east 
of the dynamic center of the LMC (see Figure 1). 



While the H I and stellar components of the LMC are 
known to extend to large distances ijKim et al.l Il999t 
Ivan der Mar el & Cioni 2001), the Ha emiss ion tends to 
be more concentrated (fCaustad et al.ll200H) . N70 IS one 
of the most outlying bright emission-line objects, located 
in an isolated region of the disk far from supergiant 
shells (IMeahiirnl l198m. H II regions, and bright diffuse 
Ha emission. The emission nebula is centered on an OB 
association (LH114 ) consisting of at least nine stars of 
type BO or earher l)Oevlll996a^ . An assumed Salpeter 
initial mass function predicts one additional ^ 60 Mq 
member which has presumably already exploded as 
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Fig. 2. — The N70 sight Unes are compared to four nearby sight lines toward hot stars. The comparison stars are located in the field and 
diffuse HII emission on the eastern edge of the LMC. Angular separations from N70 are overlaid on a POSS2/UKSTU red-band image. 
The rotational center of the LMC lies 2.6° west of N70. 



a supernova l|Oevlll996hj) . The SB is consistent with 
modest shock heating, with a shghtly enhanced [S II] /Ha 
ratio, weak X-ray emission, marginally non-thermal 
radio continuum, and a mildly supersonic expan sion 
velocity 40 km s"\ see lChu fc Kenmcuttjll988l and 
references therein). N70 displays the same filamentary 
structure seen in many m ature supernova remnants 
(SNRs) f Skelton et al] 119991. but does not show any 
of the high-velocit y emission knots commonly seen in 
those same SNRs |Chulll997j) . lOev et alJ ll2000j) model 
N70 as a density-bounded nebula with a SN shock and a 
photoionized post-shock gas arranged into the Ha-bright 
filaments seen in the nebula. From optical emission 
line ratios, they calculate that shock velocity must be 
~ 70 km s~^, which is much too low to produce O VI in 
standard radiative shock models. 

We obtained FUSE observations toward four early- 
type stars projected within N70 to search for evidence of 
hot gas. The FC/S'i? instrument covers a spectral range of 
905 — 1187 A (including the O VI doublet) at a resolution 
of ~ 15 km s~^. Further details about the instrument 
and o n-or bit performance can be found in iMoos et al.l 
l|2000() and lSahnoweTall l|2000j) . 

The N70 targets are all of spectral type 09 or earlier 
and have well-characterized continua in the spectral re- 
gion near the O VI A1032 line. The weaker O VI A1038 
line is blended with strong absorption from O I, H2, and 
C II* ISM hues and is not usually usable. The FUSE 
targets are indicated in Figure 1 and the observational 
specifics are listed in Table 1. 

We gauge the O VI contribution from the LMC halo 
in the direction of N70 by selecting a group of nearby 
comparison sight lines from archival FUSE data. Though 



FUSE has observed o ver 100 early-type stars in the LMC 
(Danforth et alJ 1200 2: Danforth 200|), we require spe- 
cific attributes for optimal comparison targets. Prox- 
imity to N70 is the first requirement since we wish to 
probe a similar sight line through the LMC halo mate- 
rial; we limit our search to observations within 60' of 
N70. Furthermore, we wish to avoid contamination from 
possible local sources of hot gas and thus disqualify tar- 
gets toward H II regions, other SBs, supergiant shells, 
and SNRs. Finally we need an unambiguous stellar con- 
tinuum around 1032A and thus rule out stars of type 
BO and later. We found three 0-type stars observed by 
FUSE that match these criteria. In addition, we included 
CAL 83 (RX J0543.6-6822), a cataclysmic variable/X-ray 
binary target located in a the low-emission field south of 
N70. This target is fainter, but the continuum is well- 
defined and acceptable for use as a comparison. Figure 2 
shows that our comparison stars are located on three 
sides of N70. Thus, any large-scale trends in O VI col- 
umn density across the region should be averaged out to 
first order. We also list the observational parameters for 
these comparison targets in Table 1. 

Each of the eight data sets was reduced with Cal- 
FUSE V2.0.5 or higher, aligned, and coadded by ex- 
posure. Our analysis is based on LiFl data since this 
channel has the highest sensitivity. The exception is 
CAL 83, for which we aligned and coadded the LiFla 
and LiF2b channels, resulting in an increased S/N for 
this faint source (at the cost of a slight decrease in spec- 
tral resolution). O VI absorption lines tend to be broad 
(6 ~ 40 km s~^) compared to the instrumental profile 
(^inst ^ 10 km s^^) so this slight decrease in resolution 
is acceptable. We converted observed heliocentric veloc- 
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TABLE 2 

H2 Measurements for BI 237 



TABLE 3 
O VI A1032 Measurements 



J 


V 


6 


log AT 


J=3 


+5 ± 5 


15 ±5 


14.7 ±0.2 




281 ± 1 


16 ±2 


15.24 ±0.10 


J=4 


~ -5 


~ 16 


< 14.2 




281 ± 1 


14 ±2 


14.90 ±0.12 



Target 



V limits^ 



log Na 



ities to the local standard of rest (LSR) with a uniform 
correction of dlsr = vq — 15.7 km s~^. The wavelength 
scale of each sight line was then calibrated such that the 
strongest component of the narrow Si II A 1020. 70 line 
was at v\sr = km s~^, an adjustment of much less than 
a 20 km s~^ resolution element in all cases. Data from 
the LiFlb segment were aligned with the LiFla segment 
by cross-correlating Fe II Al 144.94 absorption with Si II 
A1020.70. Again, this involved a shift of less than a res- 
olution element in all cases. 

Absorption-free regions of stellar continuum within a 
few Angstroms of 1032 A were used to normalize the 
spectra with low-order Legendre polynomials. The O VI 
data and continuum fits are shown in Figure 3 for both 
the N70 and comparison sight lines. The shaded regions 
denote the ±la range of continuum placements. Note 
that the CAL83 continuum is generally well-defined ex- 
cept in the region of the stellar O VI emission feature 
redward of the LMC absorption line. We choose to fit a 
smooth curve to the blue side of the emission line and 
ignore continuum longward of this feature. 

Several H2 lines can potentially contaminate the 
O VI A1032 profile, specifically absorption at Galactic 
and LMC velocities in the (6-0)P3A1031.191 and (6- 

0) R4A1032.349 transitions. Fortunately, absorption in 
J > 3 is less common t han for J < 3 toward the LMC 

1) Tiimlinson et al.]l2002() . The only sight line with sig- 
nificant J > 3 H2 absorption is BI237. We measured 
apparent column densities, line widths, and velocities in 
nine different unblended, J = 3, 4 H2 lines in the BI 237 
data and averaged the results for each J-state and veloc- 
ity (see Table 2). A model H2 spectrum was produced 
based on these values and convolved with the FUSE aper- 
ture function. The resulting model was divided from the 
observed O VI+H2 profile around 1032 A in the BI237 
data to derive a corrected profile. The H2 model is shown 
in the BI237 continuum fit in Figure 3. 

The Galactic {v < 150 km s^i) and LMC (w > 
150 km s~^) absorption components are distinct and 
well-resolved by FUSE for all sight lines even though 
there may be multiple absorption subcomponents hid- 
den in each system. We measure the apparent opti- 
cal depth (AOD) column density (iSavaee fc SembachI 

f9l') and profile-weighted mean velocity and line width 
cmbach & Savage 1992) for both Galactic and LMC 
absorption (Table 3). Saturation is a concern with the 
AOD method and the apparent column density Na is 
usually taken as a lower limit to the true column density 
N. In this case, we are reasonably certain that the O VI 
absorption is not saturated; the profiles are broad and do 
not approach zero. The thermal line width fethcrm — 17 
km s~^ for O VI at its peak coUisional ionization equi- 
librium temperature of T = 10^'*^ K, and this width is 



Galactic Absorption, AOD Method 



LH 114-7 


-55, 


154 


44 ± 


2 


74 


± 


2 


14.56 ± 0.02 


Sk-67°250 


-50, 


156 


53 ± 


3 


76 


± 


4 


14.60 ±0.04 


D 301-1005 


-66, 


153 


39 ± 


1 


68 


± 


1 


14.49 ± 0.04 


D301-NW8 


-53, 


155 


39 ± 


2 


63 


± 


4 


14.38t»:« 


BI272 


-55, 


174 


49 ± 





74 


± 


1 


14.47 ±0.03 


Sk-67°266 


-58, 


153 


33 ± 


1 


66 


± 


2 


14.41 ± 0.04 


BI 237 


-73, 


124 


31 ± 


1 


63 


± 


1 


14 47+0.06 
^^■^'-0.08 


CAL 83 


-72, 


137 


21 ± 


1 


56 


± 


5 


14 43+0.08 


(N70> 






44 ± 


7 


70 


± 


6 


14.51 ±0.10 ±0.04*= 


(conip) 






34 ± 


12 


65 


± 


7 


14.45 ± 0.03 ± 0.06'= 



LMC Absorption, AOD Method 



LH 114-7 


154, 


350 


243 ± 


2 


61 ± 2 


14.45 ± 0.02 






Sk-67°250 


156, 


340 


237 ± 


3 


59 ±3 


14.57 ±0.05 






D 301-1005 


153, 


323 


241 ± 





54 ± 1 


14.45 ± 0.04 






D301-NW8 


155, 


326 


248 ± 


1 


48 ±4 


14 45+0.05 

14.40_(J Qg 






BI272 


174, 


328 


245 ± 





45 ± 1 


14.30 ±0.03 






Sk-67°266 


153, 


325 


242 ± 


1 


41 ±6 


14 oq+OM 






BI 237 


124, 


336 


232 ± 


1 


70 ±2 


14.39 ± 0.03 






CAL 83 


137, 


287 


212 ± 





44 ±2 


14 30+0.07 






(N70> 






242 ± 


5 


56 ±6 


14.48 ± 0.06 


± 


0.04'= 


(comp) 






233 ± 


15 


50± 13 


14.29 ± 0.14 


± 


0.05'= 




LMC Absorption 


Profile Fit Method 






LH 114-7 






241 ± 


2 


61 ±4 


14.46 ± 0.03 






Sk-67°250 






234 ± 


4 


58 ±6 


14 61+° °* 






D 301-1005 






241 ± 


2 


52 ±4 


14 46+0.04 

i4.4D_Q Q3 






D301-NW8 






247 ± 


3 


49 ±4 


14 47+0.06 
J-*-^'-0.05 






BI272 






240 ± 


2 


39 ±3 


14.30 ±0.04 






Sk-67°266 






239 ± 


3 


38 ±4 


14.09 ± 0.06 






BI 237 






240 ± 


6 


72 ±9 


14.40 ±0.05 






CAL 83 






222 ± 


7 


34 ±9 


14 37+0.12 
-0.13 






(N70> 






241 ± 


5 


55 ±5 


14.50 ±0.07 


± 


0.04'= 


(comp) 






235 ± 


9 


46± 18 


14.29 ± 0.14 


± 


0.07'= 


^Units are km s ^. 
''First quoted uncertainty 


is 1(T deviation of 


mean values. 


Second 



uncertainty is mean of measured errors. 



comparable to the FUSE resolution element. 

Uncertainties in our measurements arise from contin- 
uum placement error and velocity integration limit un- 
certainty. To properly address this, we vary the contin- 
uum placement by ilcr and adjust the central velocity 
limit by ±10 km s~^. These errors are added in quadra- 
ture to give confidence limits on our AOD measurements. 
Finally, we calculate ensemble averages and standard de- 
viations for all quantities (Table 3). 

As a check on our AOD measurements, we fit the nor- 
malized spectra with a series of Voigt profile components 
convolved with the FUSE instrumental profile with free 
parameters v, b, and N. Component centroids were 
picked interactively and allowed to vary ±20 km s~^. 
Doppler width b is allowed to vary between 10 < < 100 
km s~^. Uncertainties are the quadratic sum of fit un- 
certainties and high-low continuum placement as noted 
above. The Galactic O VI absorption is best fit with 
a pair of Voigt components at ulsr 0,65 km s^^ for 
most sight lines. However, these two components are too 
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-600-400 -200 200 400 600 800 -400 -200 200 400 600 800 

vlsr (km s-i) vlsr (km s"!) 

Fig. 3. O VI continuum fits for the N70 (left) and comparison (right) sight fines. All data are from the LiFla chamicl except CAL83 
for which LiFla and LiF2b have been aligned and coaddcd. The N70 sight lines were observed through the 4" X 20" MDRS aperture while 
the comparison sight lines were all observed with the 30"x30" LWRS aperture. The absorption near zero velocity is attributed to Galactic 
OVI, and that near +210 km s^^ is due to OVI in the LMC. Continuum fits (dashed) and ztlcr uncertainties (shaded) are low-order 
Legendre polynomial fits to absorption- free regions of the spectrum around 1032 A. The BI 237 continuum includes modeled H2 absorption 
as described in the text. 



closely blended for strong constraints on their individual 
fit parameters. The sum of the two components is gener- 
ally consistent with the AOD A^ovi measurements, but 
the resulting errors from the Voigt fit are very large and 
we do not consider these further. The LMC O VI absorp- 
tion is generally well-fit with a single Voigt component 
at VhSR ~ 235 km s^^. These fits (Table 3) are entirely 
consistent with the AOD measurements in v, b, and TVqvi 
with comparable uncertainties on all quantities and pro- 
vide a valuable consistency check in our methodology. 

3. DISCUSSION 

The measurements in Table 3 show an excess of O VI 
at LMC velocities in the N70 sight lines compared 
with the field sight lines. The mean, integrated ap- 
parent O VI cohimn density for the N70 sight lines is 
(logA^N7o) = 14.48 ± 0.06 while the comparison targets 
show (log A^comp) = 14.29±0.14. The component fit mea- 
surements are nearly identical; (logA^N7o) = 14.50 ±0.07 
versus (logATcomp) = 14.29 ± 0.14. In either case, the 
N70 targets show ~ 60% more O VI column density 
{Novi = 10^"^ °^ cm-2) at LMC velocities than the com- 
parison targets, a difference larger than the dispersion 



in cither sample. The same is not true of the Galactic 
absorption; the mean integrated Galactic O VI columns 
are (log TV) = 14.51 ± 0.10 and (logA^) = 14.45 ± 0.03 for 
the N70 and field sight lines, respectively. While the N70 
sight lines show a slight excess in comparison to the field 
sight lines, it is only by a factor of 15% and is within the 
errors. 

Kinematically, we see that the N70 O VI profiles and 
comparison profiles have consistent LMC absorption cen- 
troids ((w)n7o = 242 ± 5 km s~-^ and (-y)comp = 233 ± 15 
km s~^) and line widths ((6)n7o = 56 ± 6 km s~"'^and 
('>)comp = 50 ± 13 km s"^). The N70 absorption is 
marginally broader and at higher velocity than the com- 
parison absorption, but this is within the error bars in 
both cases. 

To illustrate the differences between SB and compar- 
ison profiles graphically, we combined the normalized 
O VI data for the N70 and comparison sight lines pixel- 
by-pixel with an error-weighted average to create com- 
posite profiles. Combining the N70 sight lines makes 
sense as they are all within a few arcminutes of each other 
and probe nearly identical lines of sight. The comparison 
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Fig. 4. — Mean composite OVI profiles for N70 (solid) and the 
comparison sight lines (dotted). The two profiles show remark- 
ably similar structure in the Galactic component, but the LMC 
absorption in N70 noticeably stronger and broader. 



sight lines are separated by more than a degree on the 
sky in some cases and probe potentially different mate- 
rial, particularly at the distance of the LMC. Still, the 
scatter in comparison sight line centroid velocities is be- 
low the FUSE resolution (comparable to or smaller than 
(t{v) in the N70 sight lines) and the data were deemed 
homogenous enough to combine into a composite pro- 
file. From Figure 4, we see that the Galactic profiles 
are very similar in shape and that the N70 profile shows 
only slightly greater depth. The LMC components are 
centered at nearly the same velocity, but show decidedly 
different average profiles. We note that these composite 
profiles are used for illustrative purposes only and that 
all quoted measurements are averages of individual sight 
line data, not the combined data. 

Is the 60% excess O VI column density observed in 
the N70 profiles compared to nearby stars due to excess 
O VI within N70 or a fluctu ation between LMC halo sight 
lines? iHowk et al.l l)2002a|) measured O VI absorption in 
12 selected sight lines toward the LMC. These sight lines 
probed many different ISM morphologies, but they found 
no correlation between Nqvi and ISM morphology as 
traced by either local Ha or X-ray brightness. Similarly, 
they found no correlation between TVovi and the spectral 
type of the star observed by FUSE (usually the most 
massive star in the local OB-association). They found 
that Nqyi varies by a factor of ~ 4 between sight lines 
with a mean of (logA'^ovi) = 14.371q;2i i^'^ deviations). 
Despite the excess O VI seen in our N70 observations, the 
mean column densities for both N70 and the comparison 
sample fall within this range. 

We have re-examined the iHowk et al.l l|200 2a'^ sight 
lines, however, and reach a somewhat different con- 
clusion. Four of the stars used by Howk et al. 
are within SBs: Sk-68°80 within superbubble N144, 
Sk-70°91 (N204), Sk-71°45 (N206), and Sk-69°191 
(N154). These sight lines show (log A^ovi(SB)) = 14.49± 
0.12 while the other eight, non-SB sight lines show 
(log7Vovi(non-SB)) = 14.29 ±0.18 (see Table 4). These 
column densities and the relative ~ 60% difference be- 
tween them are remarkably consistent with our measure- 
ments of the four N70 and four comparison sight lines, 
respectively. This is strong corroborating evidence for a 
SB origin of the excess A^ovi- Furthermore, this assess- 



TABLE 4 
LMC A^ovi Summary 



Sample A/" (logATovi) (A^ovi>/10i'* Rcf.=» 



Howk et al. sample 


12 


14.37 


+0. 14 
-0.21 


2.34 


± 


0.89 


1 


Howk et al. sample'' 


11 


14.40 


+0.12 
-O.IB 


2.49 


± 


0.78 


1 


this work sample 


8 


14.39 


±0.14 


2.54 


± 


0.72 


2 


total LMC sample 


20 


14.35 


±0.17 


2.42 


± 


0.81 


1,2 


total LMC sample*" 


19 


14.38 


±0.13 


2.51 
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''Sources: 1 - IHowk et ahl l2002alV 2 - this work 

''Excluding Sk-67°05 sight line, seefHowk et al] 12002al) 
''Sight lines probe five SBs: N70, N144, N204, N206, and N154 



ment raises hope that other SBs can be studied system- 
atically for assessing O VI absorption. 

If this interpretation is correct, it has important im- 
plications for the hot gas in galactic halos. First, the 
Howk et al. mean halo A'^ovi is "contaminated" by con- 
tributions from four SBs. Removing the SB sight lines 
(and including our four comparison measurements) re- 
sults in a lower value for the LMC halo component (Ta- 
ble 4). Surprisingly, excluding SB sight lines from the 
LMC halo survey does not appreciably reduce the de- 
gree of variability of halo O VI. The entire sample of 20 
LMC sight lines (this work plus the Howk et al. sam- 
ple) shows (logA^ovi) = 14.35 ± 0.17 while the twelve 
non-SB sight lines show (logiVovi) = 14.28 ±0.16. Howk 
et al. note that one target (Sk — 67°05) has a question- 
able continuum placement and low wind velocity. We 
note that this target shows by far the lowest iVovi in the 
combined sample and is located on the extreme western 
edge of the LMC. If this point is removed, the means 
become (logA^ovi) = 14.38 ± 0.13 for the full 19-target 
sample and (logA'^ovi) = 14.31±0.11 for the 11 field sight 
lines. In both scenarios (with and without Sk— 67°05), 
the mean decreases significantly, but the dispersion re- 
mains essentially unchanged. 

The SB sample shows a systematically higher column 
density with (logA'ovi) = 14.48 ±0.09 for eight targets in 
five SBs. While this is less variation in log A^ovi, it actu- 
ally represents a similar absolute variation in linear iVovi 
(see Table 4). If the systematically higher iVovi observed 
toward SB targets represents a the sum of halo O VI 
plus a constant SB contribution, the halo contribution 
shows no less absolute variability than the non-SB tar- 
gets. Thus the sample of five SBs probed to-date in the 
LMC (eight sight lines) shows a significantly higher mean 
iVovi than do the twelve non-SB sight lines. The abso- 
lute scatter within the samples, however, is equivalent for 
SB and non-SBs. The database of stellar LMC FUSE 
observations of the earliest-type stars has grown sub- 
stantially and a detailed analysis of additional sight lines 
would undoubtedly improve the statistics on SB/non-SB 
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Fig. 5. — Comparision of composite OVI A1032 (solid) and 
Fell A1145 (dotted) profiles for the four N70 sight lines. The Fell 
profile shows three distinct LMC components in all four sight lines 
at V = 213, 254, and 295 km s"-"^ (arrows). 



O VI. 

Howk et al. (2002a) may have been too conservative 
in their conclusion regarding the lack of correlation be- 
tween A'^ovi and ISM morphology. They probably saw 
no correlation between iVovi and Ha brightness because 
they trace different processes; Ha brightness is largely 
a function of photoionization while O VI is a product 
of coUisional ionization and/or interfaces. Both H II re- 
gions and SBs are bright in Ha emission, but only the 
latter have the conditions required for O VI production. 
The lack of correlation between stellar spectral type and 
Novi is also peripheral as SB dynamics and O VI pro- 
duction are often driven by the evolved members of a 
cluster which have since exploded as SN, not the current 
association members observed with FUSE. 

We measure a perpendicular column density in the 
Galactic halo component (logA^ovi sin|&|) = 14.25 ±0.10 
for the N70 sight lines and 14.19 ± 0.03 for the compar- 
ison stars (6 ~ —33°). Both values are consistent with 
larger studies of Galactic halo O VI in extragalactic sight 
Unes (Ho wk et al.]l2002bt ISavaee et al.ll2003t) . but a fac- 
tor of three higher than that found in the first 10 kpc of 
the Galactic halo (Zsargo et al. 2003). 

Absorption in the LMC disk is traced by low-ionization 
species such as Fe II, and the strongest Fe II component 
often indicates th e systemic velocity a t a given locatio n 
in the LMC disk l|Danforth et alJl2(W : "Danfo rth ll200^^ . 
We normalized the strong Fe II A1145 absorption for all 
four N70 sight lines and created a composite profile in 
the same manner as our composite O VI profile described 
above. Figure 5 shows Fe II in comparison to O VI. All 
four Fe II profiles show three distinct absorption com- 
ponents of increasing strength at ulsr ~ 213, 254, and 
295 km s~^. We take the strongest of these components 
(vlsr ~ 295 k m s~^) to be warm disk gas at the N7 sys- 
temic velocity. iBlades. Elliott. Sz MeaburnI l)1980() mea- 
sured optical [O II] emission at various locations in N70 
and see a strong emission peak at the same velocity at 
both eastern and western edges of the nebula. We in- 
terpret the middle component {vlsr = 254 km s~^) as 
warm, photoionized gas in the approaching side of the 
SB, consistent with the quoted expansion velocity of ~ 40 
kms~^. Receding emission, presumably from the far s ide 
of N70, is seen by Bla des. Elfiott. fc Meaburn ' ("1980) at 
u ~ 311 km s~^ toward the center of the SB. We do not 
see any Fe II absorption near this velocity nor do we see 



O VI absorption above the systemic velocity which sup- 
ports our claim that the FUSE target stars are located 
within the SB, not behind it. 

The weak i^lsr ~ 213 km s~^ component is similar to 
that seen in a number of other LMC sight lin es with no 
appa. rent correlation to LMC structure (Dan forth et alJ 
120021) . We interpret this as either a high- velocity cloud 
along the line of sight, or a series of neutral LMC halo 
c omponents. 

Ho wk et al.1 ([2002a) found that LMC O VI profiles were 
systematically shifted ~ —30 km s^^ with respect to Fe II 
absorption. This is consistent with either a halo that 
corotates with the LMC disk, or a vertical outfiow of hot 
gas from the disk at v± ~ 40 km s~^. Unfortunately, the 
~ 40 km s~^ expansion velocity of N70 makes SB O VI 
absorption kinematically indistinguishable from typical 
halo O VI absorption. Comparing absolute component 
velocities to better than ~ 10 km is difficult with 
FUSE data, but we note that the O VI centroid veloc- 
ity of 242 ± 5 km s"-^ is somewhat lower than the 254 
km s^"'^ Fe II component assumed to represent the ex- 
panding shell. This is also Av ~ 60 km s~^ from the 
assumed systemic velocity. The comparison sight lines 
show VQYi = 233 ± 15 km s^^. 

There is a significant variability in A'^ovi between differ- 
ent LMC sight lines (Howk et al. 2002a). However, this 
is the first study to examine closely spaced sight lines 
through a single ISM structure in the LMC. Three of the 
four N70 sight lines show A^ovi consistent within their 
error bars in both the AOD and Voigt fitting methods. 
The exception (Sk— 67°250) shows noticeably stronger 
absorption. The uniformity of the three sight lines over 
< 87" (< 22 pc at the distance of the LMC) suggests that 
O VI is distributed evenly throughout N70 on roughly 
this scale. However, the Sk — 67°250 sight line is located 
only 16" 4 pc) from the LH 114-7 sight line and shows 
an excess O VI column of ~ 0.1 dex. This angular scale 
is similar to that of Ha filaments seen at the edges of N70 
and the column density difference may indicate that O VI 
is distributed nonuniformly on about the same scale. 

We now discuss the possible origins of the observed 
A^ovi~ioi* cm^^ excess column in N70. Four possibilities 
include photoionization, shock heating, turbulent mixing 
at the interface between the hot interior and the cold 
shell, or thermal conduction. Photoionization of 0^+ 
requires a harder spectrum (photons of > 114 eV) than 
even the hottest O stars can produce^. 

O VI is often attributed to shocks, and indeed observed 
SNR column densities are consistent with the N70 ob- 
servations (iVovi = 2.5 X 10" cm-2 toward S NR 0057- 
7226, Hoopcs et al. 200^. However, lOev et al. (200C|) 
infer Wshock ^ 70 km s""'^ based on optical emission line 
ratios across the face of N70. Although this is high for 
the expansion velocity of a SB, it is far too slow to pro- 
duce any significant amount of O VI. Furthermore, our 
sight lines intersect the SB shell nearly perpendicular to 
the presumed direction of expansion, yet we see O VI 
absorption near the systemic velocity, inconsistent with 
a shock front fast enough to produce O VI. 

Turbulent mixing layers involve a hot, low-density 

^ While OVI may be produced locally via photoionization by 
extremely hot white dwarf stars (Dupree & Raymond 1983), this 
is considered unfeasible in large-scale ISM contexts. 
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phase and a warm, high-density phase with a relative 
motion between them. Cooler material is entrained and 
heated by the hot flow in turbulence at all scales and OVI 
and related ions may be generated. Given the known ex- 
pansion of N70 and the convoluted morphology of the 
optical shell, such relative motions are plausible. How- 
ever, current simulations predict very low column den- 
sities (iVovi ^ 10^^ c m~^) per interface over physical 
scales of a fe w parsecs l|Slavin. S huU. & Bcgclman 1993; 
[Slavin 2006) . While multiple turbulent mixing layers 
may explain the observed O VI in the halo, there is not 
enough room for sufficient interfaces within a structure 
as small as N70. 

Thermal conduction models seem to provide the best 
explanation for the observed TVovi: but still have diffi- 
culty producing sufflcient column density. In this mech- 
anism, a hot, low-density phase (T ^ 10^ K) IS m con- 
tact with a cooler, high-density phase. Energy from 
the hot gas propagates into the cooler material as an 
evaporation front. When a significant amount of ma- 
terial reaches T ~ 10^ K, radiative cooling becomes 
dominant and the evaporation front stalls. A conden- 
sation fro nt then p ropaga tes b ack in to the hot side. 
iBorkowski . Balbus. & Fristr omI 1)19901) model a plane- 
parallel conduction interface and with O VI present in 
both evaporation and condensation front s with a col- 
umn d ensity A^ovi ~ few x 10^'^ cm~^. I Weaver et al.l 
l)1977|) model a wind-blown bubble with a conductive in- 
terface between the hot, shocked stellar wind and the 
cool, swept-up ISM and find ^ovi « 2.3 x 10^^ cm'^. 
ISlavin & Cox (1992) calculate A^ovi = 1-5 - 6.8 x 10" 
cm^^ for old SNRs, depending on a number of param- 
eters. More recent models of old cooling bubbles come 
close to A^ovi 10^"* cm^^ (J. D. Slavin, priv. comm.), 
but cannot quite explain the observed O VI excess in 
N70. It is also unclear if there has been sufficient time in 
the evolution of N70 for conductive interfaces to generate 
their maximum potential Nqvi- 

With the complexity of the nebular structure of N70 
as seen in Ha, we could invoke multiple interfaces to 
bring the total A^ovi up to the observed values. However, 
it seems unlikely that either all four sight lines pierce 
the same number of strong interfaces (such as conduc- 
tive interfaces), or that there is enough volume inside 
N70 for a statistically large number of weak interfaces 
(such as turbulent mixing layers). We conclude that a 
single thermal conduction interface is the most likely 
source of the observed hot gas, although some combi- 
nation of these options is also possible. Observations 
of other Li-hke ions (N V, C IV, Si IV) woul d be of 
great help in differentiating he ating scenarios ijSpitzed 
ll996llIndebetouw fc Shuli200l . although no instrument 
is currently available to make the needed observations. 
Photoionization will likely also play a role in the produc- 
tion of these lower ions, but many of the critical parame- 
ters (spectral type of the photoionizing stars, gas density, 
etc) are already well-constrained making the modelling 
relatively benign. 

As a further reality check, we estimate the total hy- 
drogen density associated with the O VI implied by our 
measurements. Assuming that the O VI is distributed in 



a shell of thickness I, then 



Npvi/l 
(O/H) foYi 



(1) 



Assuming an O VI fractional ion abundance /ovi = 0.2 
and a typical LMC H II region abundance {O/H) = 
2.34 X 10-* (Russell & Dooita 1990), we derive nn = 
0.75/Z cm^"^ pc~^. For an O VI shell occupying the 
outer 10% of the SB radius (^ '--^ 5 pc), we get nn ~ 
0.15 cm^'^. This value is nicely intermediate between 
the inferred densities of the photoionized shell (nn ~ 
10 cm~^, fOcv ct al. 2000) and the hot bubble interior 
(nn 10-3 cm-^ lOev fc Garcia-Seg'urall200^ . as one 
would expect from a thermal conduction interface. 

In addition to the O VI absorption seen toward the stel- 
lar targets, we can search for O VI emission from N70. 
The four N70 sight lines are in a crowded stellar field and 
were observed with the 4"x20" FUSE MORS aperture. 
The 30"x30"LWRS aperture is positioned 208" away (52 
pc at the distance of the LMC) from the MDRS aper- 
ture at a position angle determined by the FUSE space- 
craft roll angle at the time of each observation. LWRS 
aperture placements can potentially constrain O VI emis- 
sion and, in conjunction with column densities measured 
in absorption, determine the gas density novi-^ Unfor- 
tunately, the LWRS aperture locations for all four N70 
pointings fall just outside the bright Ha limb of the SB 
as shown by the dotted boxes in Figure 1. Hence, we 
instead use these observations to place an upper limit on 
any general diffuse O VI emission in the region near N70. 

None of the individual LWRS pointings show emission 
in cither of the O VI doublet lines at either Galactic or 
LMC velocities. Coadding the four LiFla data sets to 
increase the S/N, we have a total of 41.7 ks of expo- 
sure time. The coadded data show no evidence for O VI 
emission and we set a 2a upper limit for diffuse O VI 
emission of Fqvi < 1-8 x 10~^ erg cm~^ sr"* or 
^ovi < 9600 photons cm~^ s""'^ sr~*. This upper limit 
is not as sensitive as the 2a limit on diffuse O VI emis- 
sion in the loc al bubble of < 500 photons cm~^ s~^ sr~^ 
l)Sheltonll200^ bu t is ~ 20 time s lower than the O VI 
emission found bv lDanforth et a l. (2003) on the limb of 
SNR 0057-7226. A LWRS aperture placed within the ring 
of N70 could potentially show high-ion line emission and 
would be a definitive test for SB O VI. 

Finally, we note that the composite N70 and compar- 
ison profiles don't match perfectly in Figure 4 in the 
Galactic component. The N70 profile shows a slight 
excess in absorption at wlsr ~ 120 km s~^ while the 
profiles at w « km s^* match quite well. This ex- 
cess is nominally in the Galactic absorption compo- 
nent but could plausibly represent hot material escaping 
from N70 with a projected blueshift of ~ 150 km s^*. 
lOev fc Garcia-Segural l)2004D note that many SBs expand 
slower than predicted given the mechanical energy input 
of stellar winds and SN. They posit that hot gas may 
be escaping through holes in the photoionized shell thus 
decreasing the internal pressure. This process leads to 
galactic chimneys on larger scales and is thought to be 



^ For instance, IDanforth et alj <2003l) detected Fqvi = 4.1 x 
10~^ erg cm~^ sr~^ in a LWRS aperture placed serendipi- 
tously on the edge of SMC SNR 0057-7226 in an exposure of only 
11 ksec. 
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responsible for the general hot phase of the ISM. A hot 
component escaping at the observed velocity might be 
evidence for such an "O VI leak." If confirmed, such a 
model might increase the viability of shocks contributing 
significantly to the observed O VI excess from N70. 

4. CONCLUSIONS 

We have compared FUSE data for four sight lines 
through the ~ 105 pc diameter LMC superbubble N70 
with four nearby sight lines sampling diffuse and field 
regions in the LMC. This position on the eastern edge 
of the LMC is isolated from other large-scale ISM struc- 
tures such as supergiant shells, H II regions, supernova 
remnants, and other SBs. Thus, our sight lines are sub- 
ject to minimal sight line confusion. We see strong O VI 
absorption in all eight sight lines organized into two clear 
absorption components; one at —50 < wlsr < +150 
km arising in the Galactic halo, and another at 
+ 150 < WLSR < +350 km s"^ from LMC gas. We 
measure O VI absorption using both integrated apparent 
optical depth and profile-fitting methods; both methods 
yield consistent results. The Galactic column densities, 
line widths, and profile shapes are equivalent for the N70 
and comparison sight lines and are consistent with typi- 
cal Galactic iVovi seen in other surveys. The N70 sight 
lines show a ~ 0.2 dex (~ 60%) higher column density 
at LMC velocities than the comparison sight lines, larger 
than the scatter in either sample. This corresponds to 
NoYi « 10" cm-2. Furthermore, the N70 O VI profiles 
show slightly broader absorption than do the comparison 
profiles. 

Our measured iVovi values for both N70 and com- 
parison sight lines are consistent with those found by 
|Howk et al. ( 2002a) in an earher survey of LMC O VI 
absorption along twelve sight lines. While Howk et al. 
claim no correlation between A^ovi and ISM morphol- 
ogy (as measured by either stellar spectral type or Ha 
or soft X-ray surface brightness), we find that our N70 
Novi is consistent with that measured in four SB sight 
lines in their sample. Similarly, our four comparison mea- 
surements are entirely consistent with the eight non-SB 
measurements in the Howk et al. sample. We show that 
the SB and non-SB samples show distinct mean Nqvi 
but that the scatter (as expected from halo variation) is 
equivalent for both samples. This is strong evidence that 
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SBs (including N70) produce a local reservoir of O VI dis- 
tinct and distinguishable from the general, overlying halo 
absorption seen in nearly all LMC sight lines. 

We compare O VI absorption to Fe II, a reliable tracer 
of the warm LMC disk gas. The strongest Fe II absorp- 
tion traces the systemic velocity of LMC disk gas near 
N70 at fLSR = 295 km s~^ and the observed O VI ab- 
sorption is consistent with an expansion velocity of ~ 60 
km s~^. However, it is also consistent with the general 
trend of blueshiftcd halo O VI absorption with respect 
to Fe II seen bv H owk et al.. (,2002a,)- 

We consider four different ionization mechanisms (pho- 
toionization, shocks, turbulent mixing layers, and con- 
ductive interfaces) and conclude that thermal conduction 
between the interior hot. X-ray producing gas and the 
cool, photoionized shell of N70 is the most feasible pro- 
duction mechanism for the observed O VI. The inferr ed 
shock velocity (-Wshock ~ 70 km s~^. lOev et al.ll2000() is 
insufficient to ionize the SB shell to the observed level. 
We calculate the total hydrogen density implied by our 
iVovi measurements and find nn 0.15 cm^"^ for a 5 pc 
O VI shell, consistent with expected values. 

Our LWRS apertures lie just outside the shell of N70 
and show no O VI emission to a fairly sensitive level 
(^bvi < 1-8 X 10"'^ erg cm~^ sr~^). Future detec- 
tions of diffuse O VI emission toward N70 will conclu- 
sively demonstrate the existence of O VI within the SB. 
In addition, emission-absorption measurements can con- 
strain the density and distribution of hot gas within N70. 
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